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ABSTRACT 



The decision problem of the allocation of available weight and 
space in an air-dropped anti-submarine torpedo to fuel^ explosive^, and 
to attack speed capability above some preset minimum is studied for the 
case in which alternatives under consideration do not differ with respect 
to their effect on weapon reliability and to their effect on enemy 
countermeasures capability * Analytical models are developed which 
relate (l) probability of acquisition to endurance for the case of 
circling search at constant depths (2) probability of hit given acqui- 
sition to endurance and attack speedy (3) probability of hit on a 
close-in attack run to attack speedy and (k) probability of kill damage 
given a hit to amount of explosive usedo 
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CHAPTER I 



INTRODUCTION 

The Problem 

An air-dropped anti-submarine torpedo may be characterized by a 
long list of parameters o Specification of parameter values in a 
specific weapon is made such that the weapon may achieve a maximum 
of effectiveness in an actual tactical situation in wartime 0 The 
state of the art in weapon design will determine the values of many 
of the weapon parameters* For example* the transducer parameters 
may be prespecified by the status of some continuing torpedo sonar 
development program. Cost considerations may also be of great 
importance in the evaluation of alternative proposals for the 
design parameter values. For example* a choice between two propul- 
sion systems may hinge on a difference in development costs necessary 
to produce a reliable weapon. However* for torpedoes and more 
especially for air-dropped torpedoes* the restrictions of weight 
and of physical dimensions (or volume) resulting from the varied 
mix of capabilities of the prospective weapon delivery systems are 
also very significant in the weapon design. 

Consider the restraints of weight and volume. One of the design 
parameters which has a great effect on weapon weight and volume is 
the maximum operating depth. An increase of a few hundred feet in 
the operating depth requirement will require a more substantial 
structure. This in turn may result in a decreased weight and volume 
to be made available for other weapon components. However* intelli- 
gence reports or estimates of the depth capabilities of prospective 



submarine targets in the time period for which the weapon is to be 
designed may prompt the decision maker to fix the design operating 
depth as an input to the design problem e Similarly the advance speei 
fication of a minimum attack speed capability may be prompted by 
intelligence estimates and will certainly act to place an upper limit 
on weight and volume available for other weapon components*, Suppose 
that a study of specifications from higher authority^ engineering 
capabilities to produce the various weapon components^ cost considera 
tions or some appropriate optimization procedures have resulted in 
specification of values for all weapon design parameters except 
1* Amount of fuel 

2. Amount of explosive 

3. Maximum speed (above some preset minimum value ) c 

Further, suppose the restraints on weight and volume have been deter- 
mined such that if the amount of fuel and explosive and the maximum 
speed capability above some specified minimum are considered as a 
weapon payload, then the available payload weight and volume are 
provided to the decision maker 0 The decision problem of allocation 
of available payload weight and volume to amount of fuel p amount of 
explosive, and to maximum speed capability above some specified 
minimum is the problem considered in this paper 0 



2 



The Approach 

An initial approach to the allocation problem is to set certain 
minimum limits on the parameter values . The minimum value of maximum 
(attack) speed has been specified 0 There must be sufficient fuel to 
provide some minimum probability of target acquisition 0 (in the case 
of continuous circular search at a constant preset depth as is con- 
sidered in this paper a minimum probability of acquisition might 
result from a dive to search depth and completion of one complete 
search turn.) Then there must be sufficient fuel for at least one 
attack under some assumed conditions of acquisition range and target 
course and speed. Finally there must be at least enough explosive 
to result in kill damage if a hit occurs at the most favorable depth 
and point of contact on the target surface. Having determined these 
allocation minimums there remains the problem of allocating the 
balance of available payload weight and volume. 

A second step may determine a maximum value of attack speed. 

This maximum may be determined from a study including the probable 
operating depths the characteristics of the guidance system^, and the 
generated self-noise level and its effect on the ability of the sonar 
to maintain contact with the target on the attack run. 

What is the next step? How does the decision maker decide between 
enough additional fuel for one reattack and enough additional explosive 
to produce kill damage over a given additional element of target 
surface? What measure of effectiveness is used? 

This paper proposes a procedure which may provide the decision 
maker with a quantitative analysis or comparison of the tradeoffs 
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involved in alternatives of allocation among fuel, explosive, and 
additional attack speed capability 0 The tradeoffs are measured in 
terms of probability of kill as a measure of effectiveness * Probabil- 
ity of kill is dependent upon weapon reliability and upon the 
capability of the enemy to employ countermeasures (in addition to 
maneuvering)* It is assumed in this paper that weapon reliability 
and enemy capability to employ countermeasures (in addition to maneu= 
vering) are both independent of the amount of fuel and explosive and 
of the weapon attack speed within the range of values between the 
preset minimum and some determinable maximum* This enables removal 
of reliability and enemy countermeasures capability (less maneuvering) 
from the problem and probability of kill may then be factored into 
three independent probabilities j those of probability of acquisition, 
probability of hit given acquisition, and probability of kill given 
hit* An optimization is made of probability of kill as a function 
of the fuel, explosive, and additional attack speed capability with 
weight and volume restraints* This does not provide a number for 
use as an operational value of kill probability* However the 
dynamic programming methods used provide a solution which specifies 
the optimal allocation of weight and volume* 

It is emphasized that a quantitative comparison of the tradeoffs 
involved in the allocation process may be obtained only by some 
common analysis of the three study variables as they act upon the 
events of acquisition, hit, and kill* 
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The dynamic programming optimization model developed generates 
an input requirement for (l) probability of acquisition as a function 
of available fuel^ (2) probability of hit as a function of available 
fuel and attack speed and (3) probability of kill damage as a function 
of explosive weight. These requirements may be satisfied in a number 
of ways including 

1. analytical derivations in terms of obtainable or estimable 
inputs 

2. simulation (war gaming) 

3 0 extrapolation of empirical data from operational or 
development tests 
4. a combination of the above 

The model user may choose any of these methods so long as the 
resulting inputs are compatible with the allocation model . However^ 
any of the methods used will require an assumption about the tactical 
situation with respect to which the optimization is to be made c 
Consider the actual tactical situation in which the weapon will be 
used. Assume that methods exist to determine the optimum values 
of all weapon parameters^ given that the weapon will be used in a 
completely specified tactical situation. The optimum parameter 
values determined will be different for each different tactical 
situation considered. The space of different tactical situations 
is infinite. However a single set of parameter values must be spe° 
cified for the weapon^ ultimately to be produced in quantity. By 
what criteria is the specific tactical situation to be selected for 
use in the optimization of the weapon parameters? The probability 
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that the weapon will actually be used in the tactical situation selected 
is zero. The following two criteria will be used in this paper in the 
analytical derivation of the three required probabilities M 
Criterion 1; Select an unfavorable tactical situation which is not 
extreme. This may be thought of as a modified f maximin l criterion . 

A maximin criterion implies optimization of the weapon parameters to 
maximize effectiveness for that tactical situation which offers 
minimum probability of weapon success. The extreme situation is 
modified in this paper because the low probability of kill resulting 
from assumption of an unfavorable situation in the extreme would 
undoubtedly distort the relationships between the study variables which 
would apply in a greater percentage of cases. In general a weapon 
which is optimized for an unfavorable situation will be able to cope 
with a more favorable situation for which its parameter values were 
not optimized. 

Criterion 2 : Select a tactical situation which is sufficiently 

simplified as to be within the scope of this paper. Progress in 
work of this kind is more often made by starting with something 
relatively simple and adding sophistication after experience is 
gained with the model. 
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Relation of this work to previous work 

The allocation problem has been studied before by R, C e Brumfield 
in references 1 and 2 in which algebraic expressions were derived 
relating total weight and volume of a weapon of arbitrary speed and 
range to performance characteristics anti component characteristics. 
This enabled determination of optimum running speed and range for 
stern chase vehicles under conditions of minimum propulsion system 
weight, minimum total weight of propulsion system and associated 
shell, and minimum propulsion system volume 0 The approach in this 
paper is essentially different from that of references 1 and 2 
although the problem is much the same* The previous work optimized 
total weight and volume and attack speed for the event of a single 
attack run, given explosive weight, a length of run in the acquisition 
phase, and acquisition range. This paper optimizes an allocation of 
payload weight and volume to fuel, explosive, and attack speed for 
the events of acquisition, hit (attack and reattacks) and kill 
given the maximum allowable payload. 

A great deal of work has been done in the study of various 
phases of the torpedo problem. References 3 and k form the background 
on which the development in Chapter III is based. Reference k 
develops a continuous volume scan rate for a single circular search 
turn. In this paper, in the first three sections of Chapter III, 
a volume scan for a single ping on the first search turn is developed. 
This is a modification of the continuous scan rate to the discrete 
case where the discrete interval is a single ping cycle. However, the 
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development in this paper is continued to obtain a probability of 
acquisition not only for a single circular search turn but for con° 
tinuous turns until fuel exhaustion * This further development is 
made possible because a tactical situation is postulated in accordance 
with criteria one and two e The postulation is Chapter III limits 
the usefulness of the specifics of the development to the search 
pattern considered but is required in order to obtain the probability 
of acquisition input to the allocation model 0 The logic of the 
development may be modified as appropriate to provide for other 
postulated tactical situations in accordance with requirements to 
study a particular search pattern 0 The study of search turns subse~ 
quent to the first and until fuel exhaustion is necessary to provide 
probability of acquisition as a function of fuel weight throughout 
the range of the argument * 

The development in Chapter IV of probability of a hit in any given 
pass at the target has no known precedent in the literature and is 
presented as an alternative way of looking at the hit probability 
problem which may be of interest « 

It is emphasized that Chapters III, IV, and V are samples of 
possible methods of attacking the problem of finding inputs to the 
allocation model of Chapter I which are specifically designed to 
permit examination of the tradeoffs between alternative allocations 
of payload weight and volume among the three study variables 0 
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Structure of the paper 

Chapter I contains the development of a mathematical model of 
the torpedo problem in a form such that the dynamic programming 
techniques of Chapter II of reference 5 ma y be applied to solve the 
allocation decision problem. Chapters III^ IV^ and V contain the 
development of the three required inputs to the allocation model of 
Chapter 11 % probability of acquisition^, probability of hit given 
acquisition and probability of kill given hit. Appendix A is a 
description of the computer flow chart notation to be used in the 
flow charts of the other appendices. Appendix B is the computer 
flow chart of the dynamic programming allocation model of Chapter I 0 
Appendices C 1 through C 8 contain the mathematical details of the 
developments in sections I through VIII of Chapter III plus the 
corresponding computer flow charts. This mathematical supplement 
to Chapter in is included in appendices to enable emphasis of the 
logic flow in Chapter III with a minimum of mathematical detail. 
Appendix D contains the computer flow chart for Chapter IV. 
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Summary 

A mathematical model is provided which enables optimization of 
the decision problem of the allocation of available weight and space 
in an air-dropped anti-submarine torpedo among three study variables® 
amount of fuel* amount of explosive and attack speed capability above 
some preset minimum. The model is restricted to the study of alterna- 
tive allocations which do not differ with respect to their effect on 
weapon reliability and to their effect on enemy countermeasures 
capability (other than maneuvering capability.) 

Mathematical models are developed which relate (l) probability 
of acquisition to endurance for the case of circling search at 
constant depth* (2) probability of hit given acquisition to endurance 
and attack speed* ( 3 ) probability of hit on a close-in attack run to 
attack speed* and (U) probability of kill damage given a hit to amount 
of explosive used. These four models are developed as specific 
examples of inputs needed to be compatible with the requirements of the 
allocation model. Computer flow charts are provided for each model 
except the last to facilitate programming in any desired computer 
language and also as an alternative method of study of the development 
of the model logic. 
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CHAPTER II 



THE ALLOCATION MODEL 

Consider the design of an air-dropped anti-submarine torpedo such 
that the amounts of fuel and explosive to be included and the capabi- 
lity of attack speed in excess of some predetermined minimum value 
may be considered together as a payload for which there is available 
a predetermined maximum of weight and volume 0 Consider the employ- 
ment of the torpedo from the time it splashes into the water until 
the time it explodes in the proximity of a submarine target as an 
ordered sequence of the following three events each of which must 
occur in succession for accomplishment of the torpedo missions 

1. Acquisition after a search run following some preset pattern 

2. Hit after an attack (and if necessary a series of reattacks) 

3. Kill damage resulting from explosion ot the weapon in the 
proximity ot the target 

Divide the study variable, amount ot fuel, into two separate and 
independent study variables: 

1. Amount ot fuel to be used in tne event of acquisition (event l) 

2. Amount of fuel to be used in the event of hit (event 2) 

Restrict the study variable, attack speed, to a range of values 
such that 

1. The lower limit of the range is some preset minimum value which 
is greater than the given value of weapon search speed 0 

2 . As the attack speed is allowed to vary over its range, the fuel 
expenditure rate at attack speed and at any depth under considera- 
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tion may also vary. However as the propulsion system is modified 
to allow for changes in the value of its maximum (attack) speed 
capability, the fuel expenditure rate at search speed and depth 
remains approximately constant. 

Make the following assumptions? 

Assumption 1: Weapon reliability and target capability to employ 

countermeasures other than maneuvering are both independent of changes 
(within the ranges to be considered) in the values assigned to the 
four study variables: 

1. Amount of fuel for event 1 

2. Amount of fuel for event 2 

3. Amount of explosive 

4. Weapon attack speed 

Assumption 2: The weapon is 100$ reliable. 

Assumption 3: The enemy target submarine has no countermeasures 

capability except that of maneuvering. 

Note: Assumptions 3 and 4 are made in order to justify the relation- 

ship between the study variables and the measure of effectiveness to 
be introduced later. The fact that they are intuitively incorrect 
does not detract from the validity of the optimization in the model 
because of assumption 1. 

Assumption 4: The three events of acquisition, hit and kill are 

independent events in a probabilistic sense. That is, if 
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P [xj = Probability that event x occurs 
= P (target acquisition (event 1 )] 

Pjj = P [target is hit (event 2) given that event 1 has occurred j 
Pp = P (kill damage is inflicted on the target by explosion of the 
weapon warhead (event 3) given that event 2 has occurred 



If P R = P [target kill] 

then P = P • P • P 
nen r k “ r A r H E 

Define : 

W = available payload weight (lb) 

V = available payload volume (ft^) 

= fuel weight for acquisition, event 1 (lb) 
Wg = fuel weight for hit, event 2 (lb) 



w^'s weight for modifications to the propulsion system to provide 
a maximum (attack) speed capability in excess of some preset 
minimum value (lb) 

= preset minimum value of attack speed (yd/sec) 
w^ = explosive weight (lb) 

o 

v^ = fuel volume for acquisition, event 1 (ft 3 ) 
v^ = fuel volume for hit, event 2 (ft^) 

v^ 1 = volume for modifications to the propulsion system to provide 
a maximum (attack) speed capability in excess of S q (ft J ) 

O 

v^ = explosive volume (ft J ) 

w^ = payload weight for hit, event 2 (lb) 

i it 

w 2 = w 2 + w 2 
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v 2 = P a yload volume for hit, event 2 (ft*) 

V 2 ‘ V 2 + v 2 

A w = increment of weight of the any of the weight variables^ all of 
which are varied in increments of equal size (lb) 
gi( w i) = the relationship between the i th weight variable and 

I V 0 

the i th volume variable^ i « 1, 2, 2 , 3 

t t t 

g*^) = v £ ; i » 1, 2, 2 , 3j g L (w) = g^(w) 

f 1 C w i ) = natural logarithm (in) of P^ expressed as a function of 

W V ) = In P„ expressed as a function of w* and w“ 1 
c. c. d H d. 2 

f,(w ) = In P expressed as a function of w, 
j J r 3 

3 

F = In P = XI f i 
i = 1 1 



I _i 

w , w = lower and upper limits respectively on the allowable range 
of fuel weight (lb) 

i i i t 

W > W = lower and upper limits respectively on the allowable range 
of weight for modifications to the propulsion system to 
provide a maximum speed capability in excess of (lb) 



* i i til 

W W 

— 9 = lower and upper limits respectively on the allowable 

range of explosive weight (lb) 

The concept ot the lower and upper limits on the ranges of the 
study variables is as follows* It will be obvious to the design 
engineer that very small values or either fuel or the explosive 
weights will result in a negligible probability of success* Thus 
there is no need to examine the vicinity of the extreme values of 
zero and W . However the values of the limits may be chosen as 
close to zero and W as is desired* In the case of W ° 0 9 the value 
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chosen may well be zero as this implies that the considered minimum of 

8 

attack speed is to be S q , the preset minimum. In the case of W the 
value chosen will probably be well under W as the maximum of the range 
of values of attack speed is limited by considerations of resultant 
guidance problems, excessive self-noise levels, and so forth . The 
hit event, event 2, contains two study variables which must be con- 
sidered together. For events 1 and 3 only one study variable is con- 
sidered per event (or programming stage). Thus for each value of 
fuel weight (w^) there is associated a single value of acquisition 
probability P (w 1 ) and a single value of volume v, * g^Wp). 

Similarly for values of explosive weight w^« 

However 

W 2 = W 2 + W 2' 

T 2 = T 2 * t 2 * = 8 2 ( w 2 ) * 8 2 * (w 2 5 

This requires that for each value of weight (w^) there be an 

associated set of values of hit probability ^( W 2 ^ an< ^ volume 

v . one value for each combination of values of w' and w'* which 
2 7 d d 

♦ tt j ® 8 

sum to the considered value of w^. A suboptimization of w^ and w^ 
to provide a single optimal value of ) for each w^ as an input 

to the allocation model would neglect the associated volume 
restraints. Therefore index the weight w^ by k(k - 1, K) over 
its range 

min £w' , W"”^ ~ v 2 “ mijl ^ W ‘ + ’ » w | 
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and define 



w k = k the value of weight to be distributed between w^ and 



w“; k = I,-.-, K (lb) 



Ju (lb) 



W 2jk = •) bb value of w^ associated with >5^2 j - 1^, ” <, 

W 2jk = J value of w^' associated with w^j j s (lb) 

P k = J th value of associated with W^j j s 1 #’°% 



V 2jk = J bh value of v^ associated with w^j j - 1* J k ) 

k = 1, •••, K 



V 2jk = -i bb e value of v” associated with w^j j - 1 



> "*> * 1^1 



v = v* . + v'\, = j the value of v. associated with 
2 jk 2 jk 2 jk J 2 



w 2k i J = 1> \ ( ft ) 



The available payload W and the upper and lower weight limits or 
the weight increment A w are so chosen that all the weights are 
divisible by Aw. Then there will be 



W* W» 
Aw ” Aw 



considered values of w^ denoted w^j r = 1, 
there will be 



min { W» + W" , W J- min £ W» , W« *j- 



Aw 



Aw 



considered values of w^ denoted w^J k = 1 



- K 

K 
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and there will be 



W" W« " 



Aw 



Aw 



= S 



considered values of denoted w_ s s = l s * * ° 9 S 

3 3s" * 

The functions g^w^, g^Cw^), g^' (w^' ), g^w^)* f^Wj), and 
^3^ W 3^ Sre assumec * to non ~d ecrea sing . That is* an increase in 
the argument will not decrease the value of the function 0 In the 
case of fg(wgj w^ 1 ) it is assumed that any increase in either of 
the arguments while maintaining the other fixed will not decrease 
the value of the function. 

The problem is to maximize F = In 
subject to 

< W 

V 1 + V 2 + V 3 — V 

Consider the new problem to maximize 

p» = F - X j^g 1 (w 1 ) + g£(w£) + g£’ (w£>) + g 3 (w 3 )] 



subject to 



+ w 2 + <W 






g i^ w i^ 1 = 2 » 2 9 3 

where X is a non-negative real number (Lagrange-multiplier ) , 



Define the recurrence relations! 



h 1 (w) = f x (W) - X g 2 (W)' if W = ^ T = 1, •••, 



R 



= - o o 



otherwise 
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h_(W) = max 

^ w • k = 1 

w 2k» k ±, 



... K (W V 'V ' XT 2Jk] 

C + h l (w “ w 2 k ) 



where 



f 2^ u 2Jk’ u 2j]2 111 P Hjk “ F H , "2jk I U 2jk^ 

T 2jk = * *2’ (w 2Jk> 



3 (W) = . 3 Tl, .... s [W * WV * h 2 (W -%>| 



The dynamic programming methods of Chapter II of reference 5 are 
next applied to the recurrence relations h^, h^ and h^ for some fixed 
X to obtain a solution for the optimum value of h^(w) “ F* and for the 
corresponding optimum values of w ^ 8 9 and 

Next the corresponding optimum values of v^ and v^ are com- 

puted using the computer determined optimum values of 9 and < 

The sum (v^ + v^ 4* v^) is then compared with the available payload 
volume V. Probably v^ + v^ + v^ f V. 
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Then if 



v ! + v 2 + v 3 > V 

the value of X must be increased and the dynamic program repeated 
until the minimum value of X is found which results in 

V 1 + v 2 + v 3 <V 
If v^ + v 2 + v^ <V 

the value of X is decreased until the minimum value of X is found 
which results in 

V 1 + v 2 + v 3 <V 

The optimum value of h^w), w^, w^' and w 3 determined for 

such a X specify the solution values of the three study variables 0 
The Lagrange-multiplier, X $ may be thought of as a cost 
coefficient . If: X = ° 

v l + v 2 + V 3 ~ V 



then the volume restraint may be thought of as having cost nothing 0 
That is, the measure of effectiveness has been optimized solely with 
respect to the weight restraint and the result has met the volume 
requirements* On the other hand, if X - o 

’l*"!* v 3 > v 

then some positive value must be assigned to X which will result 
in meeting the problem volume restraint Vj^ + + v^^V. 

Ultimately the subsequent process of optimization with some 
positive value of X will result in smaller optimum values of one 
or more of the weight variables with a correspondingly reduced value 
of the measure of effectiveness. 



19 



This reduction in the solution value of the measure of effective- 



ness may be considered as the cost of having a volume restraint 0 Of 
course if weights and volumes are interchanged in the model the 
resulting X will be a measure of the cost of the weight restraint 0 
If the solution value of probability of kill is desired then 



The computer flow chart to accomplish this dynamic programming 
procedure is contained in Appendix B« 

It may be of interest to allow the values of the available 
payload weight and volume (w and V ) to vary in the vicinity of their 
preset values. Then the sensitivity of the effectiveness measure ^ 
to such small variations may be examined in order to determine 
the possible value of modifications to the weapon delivery systems to 
allow for a different (more effective) load capability „ 

Input Requirements for the allocation model 
W => Available payload weight (lb) 

V = Available payload volume (ft^) 



W , W = lower and upper limits respectively on the allowable range 
of fuel weight (lb) 

t i t 

W , W = lower and upper limits respectively on the allowable 



P k = exp 




h^(w) + X (v L + v 2 + v 
, v^, v^ and v^ are solution values. 




where h^(w) 



3 



range of weight for modifications to the propulsion 
system to provide a maximum (attack) speed capability in 
excess of a minimum preset value (lb) 
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W , W = lower and upper limits respectively on the allowable 

range of explosive weight (lb) 

t v 

gl( w l); i = 1, 2 , 3 = the relationships between the i th weight 

variable and the i th volume variable; 
i v 

i = 1, 2 , 3 

Pa( W i) = probability of target acquisition expressed as a function 
of the fuel weight for acquisition 

P„(w * ) w* 1 ) = probability of hitting the target given that target 

n 2 2 

acquisition has occurred expressed as a function of the 
fuel weight for hit and of the weight for modifications 
to the propulsion system to provide a maximum (attack) 
speed capability in excess of S q 

P (w ) a probability of target kill given that a hit has occurred 

E J 

expressed as a function of the explosive weight 
S q = minimum preset value of attack speed (yd/sec) 

£w =s size of the increments of weight of the weight variables^ 
all of which are varied in increments of equal size (lb) 

The payload weight and volume are determined by the capabilities 
of the weapon delivery systems. The limits of the weight variables 
are determined by a preliminary study which reflects limiting design 
considerations and reduces the possible ranges of the weight 
variables to reduce computer time. The weight~volume relationships 
are determined by engineering considerations reflecting the state 
of the art. The minimum attack speed is determined by intelligence 
estimates or CNO requirements. The size of the weight increments 
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is determined by considerations of required accuracy of results and 
available computer time. The probabilities of acquisition, hit and 
kill may be obtained in a number of ways including 

1. analytical derivations in terms of obtainable or estimable inputs 

2. simulation (war gaming) 

3. extrapolation of empirical data from operational or development 
tests 

4* a combination of the above 

The remaining chapters provide samples of possible analytical methods 
of providing the three required probability inputs which are specifi- 
cally designed to be compatible with the allocation model requirements 
It is possible to extend the model to provide for cost restraints 
The procedure would require a second Lagrange-multiplier in the 
function to be maximized (f 1 )<> Additional inputs required would be 
the relationsuips between the costs and the corresponding weight 
variables. It is considered, however, that cost restraints are not 
significant to the problem considered in this paper. Variations in 
the amount of fuel and amount of explosive have little effect on the 
costs involved. Even changes in attack speed within the limits 
appropriate for consideration will probably have negligible cost 
effect when compared to the total cost of developing, producing and 
operating a weapon in the course of its life cycle. 
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CHAPTER III 



PROBABILITY OF ACQUISITION 

Consider a tactical situation in which a weapon is delivered by 
some air delivery weapons system to a splash point on the surface of 
the water in the vicinity of a submarine target and make the following 
assumptions . 

Assumption 1; The splash point defines the center of the top of a 
right circular cylinder called the target volume cylinder. The 
probability is quite small and is assumed zero that the target is out= 
side of the target volume cylinder. The target is placed at random 
within the target volume cylinder. 

Assumption 2: The weapon dives to its assigned search depth on a 

straight line track and commences a right circular search turn at 
constant depth, speed, and turn rate immediately upon reaching 
search depth. The weapon searches in the active acoustic mode 
commencing its first ping upon commencement of the search turn. 

The weapon continues its turn until target acquisition or fuel ex= 
haustion. 

Assumption 3: The weapon does not acquire a false target. 

Assumption k: The search turn of the weapon is approximated by 

rotating the weapon about a vertical axis through the point at which 
the weapon commences search. 

Assumption 5: While the weapon is diving to search depth, the radius 

of the target volume cylinder increases at a rate equal to the target 

cruise speed. 
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commencement of the weapon search 



Assumption 6: The target is alerted on 
turn by the first 'ping' of the weapon. Thereafter the target opens at 
constant depth and speed radially from the projection of weapon position 
on the horizontal plane through the target. 

Assumption 7; During search the radius of the target volume cylinder 
increases at a constant rate equal to the target speed. 

Assumption 8: The 'definite range* law as developed in Chapter II of 

reference 6 applies to acquisition. Thus on a given radial during a 
ping cycle under specified conditions there is a range beyond which 
acquisition probability is zero and within which acquisition probability 
is one. Regions closer than the 'definite range' limits are considered 

scanned „ 

Assumption 9: During search the weapov makes a series of independent 

glimpses for the target, one for each ping cycle. 

Assumption 10: Target cylinder volume is considered constant during 

each ping cycle at the value valid at the initiation of the ping. 
Assumption 11: The water medium is considered to be isotropic. 

That is, all environmental parameters that affect acoustic performance 
have a constant value throughout the considered volume. 

Assumption 13: Consider planes containing the transducer and angles 

measured at the transducer. Consider the vertical plane containing the 
weapon longitudinal axis as the vertical reference plane. As the 
weapon sweeps through an angle <P in the horizontal plane it sweeps 
through a lesser angle <p » in a plane perpendicular to the vertical 
reference plane but inclined to the horizontal. Also if the inclination 
of the plane inclined to the horizontal as measured in the vertical 
reference plane is « then the inclination of the plane inclined to the 
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horizontal as measured in a vertical plane other than the reference 
plane is a lesser angle « ' . In this paper it is assumed that 

©C * = ©C 

In practical situations likely to be considered, this assumption 
results in much smaller errors than those created by treating the 
range problem in discrete steps as is done in this chapter. 

Define : 

P = p (i) = P (” acquisition in I ping cycles ] 
where P f xj = probability of event X happening 

g. = P [" acquisition during the i th ping cycle] j i - 1, •••, I 

w = weapon fuel expenditure rate at search depth and speed (Ib/sec) 

Prr = weapon active sonar pulse repetition rate for search 
s 

(pings/sec) 

w = fuel used by weapon to reach search depth, (lb) 

OSD 

w = fuel available for acquisition, (lb) 

Then, as developed in Chapter 2 of reference 6, 

(1-gj) 



P A =1 
and since 



i=l 



( W 1 “ w OSD' ) PRR s 



v 



It 



= I 



and X~ = greatest integer —X and where I 

is the number of complete ping cycles available for the event of acquisi^ 
tion, it may be said that P is a function of w L and a table may be 
computed of P A vs w x where w L is allowed to vary between W and W . 
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Of course if I is less than 1 then P A is defined to be zero* The 
probability of acquisition in a specified glimpse (g^) con ^ 

sidered as a coverage factor* The volume containing the target with 
positive probability is a region to be covered c A portion of the 
region to be covered is scanned (covered) by the weapon in accordance 
with the definite range law* Because the target has been randomly 
placed in the target volume cylinder (region of target uncertainty 
to be covered) then 

volume scanned on i th ping cycle 
®i ~ target volume during i the ping cycle 

Thus the problem is reduced to finding what portion of the volume 
of target uncertainty is scanned by the weapon on any ping* This will 
be different for every ping. 

The following development for is divided into eight consecutive 
steps * The development is amplified as appropriate and complete flow 
charts for a computer program are provided in appendices C 1 through C 8* 
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1. Available Echo Excess (B^ k ) 

Consider a vertical plane containing the longitudinal axis of the 
weapon at the time of ping transmission Consider , as in Figure 1, a 



second plane perpendicular to the first and containing the weapon 
transducer such that it makes an angle j with the weapon longitudinal 
axis, A set of such planes is considered such that the angle j is 
associated with the j th plane. In the vertical plane angles are 
measured positive counterclockwise from the reference axis. 

Design or operational considerations will make it necessary to 
prohibit reception of signals from off the transducer axis in the 
vertical plane in excess of some limting angle J, not greater than 
90 degrees. Allowing one degree between each of the considered planes j^, 

the index j will assume the range of values -J, =J + 1# ° ° °i> 

J. Thus 2J+1 planes are to be considered, one containing the trans~ 

ducer axis and J planes on either side* 



Consider increments of target range of 100 yards in the plane j. 

For an assumed target positioned in the j th plane at a range of 100k 
yards on the transducer azimuth axis projection in the plane j (point P 
in Figure l), there exists an echo excess (Bj k ) by which the signal level 



< S V 



of the echo reflected from the target exceeds the background 



masking level at the weapon transducer (M(BG) k ). That iss 



B = SLj k - M(BG) k ; j= -J> °°°» 

Assuming a definite range law of acquisition, a negative echo excess 
yields zero probability of acquisition and a positive echo excess yields 
a probability of acquisition of one. The maximum range of acquisition on 
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X direction. Line ox is the transducer axis extension. 



ox, oP, and oy are in the vertical plane. Lines oP and oz are 
in the plane j . The polar coordinates of the point P in the 
vertical plane using the x axis as a reference are (lOOkp j). 
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the azimuth axis in the j th plane is thus the range of 1 yards 

which B equals zero. At ranges less than 100K. yards there is an 
jk J 

echo excess of B above that required for acquisition on the azimuth 

jk 

axis. This excess is available for use off the aximuth axis d 

The signal level varies with assumed target range and with the 
elevation angle of the target above or below the tranducer axis in 
accordance with 

SL Jk = WS + TS - 2 H fc + 10 log ^(J) b 2 ^] 

where : 

WS = on- axis source level (db at 1 yard from transducer 
TS = target strength of target (db) 

H = one way transmission loss at a range of 100k yards (db) 
k 

t» i (j) = normalized vertical pattern function for transmitting 
pattern (i.e., the ratio of intensity at the angle j 
off the vertical axis but on the azimuth axis to intensity 
on-axis in both the vertical and aximuth planes, 
b (j) = normalized vertical pattern function for receiving pattern. 
There are several different kinds of masking noises and it is their 
total effect which is considered as the interference background masking 
level. Included are self-noise, target noise, ambient noise, and 

volume and boundary reverberation. 

There is also an enabling level designed into the weapon electrical 
circuitry to protect the acoustic system from excessive interference 
and to avoid interpretation of the interference as a signal. As used 
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here, the weapon enabling level (WEL^) includes the weapon recognition 
differential. It is a function of range (or time) and decreases with 
increasing range. The background masking level (m(BG ) r ) is the larger 
of the two: interference background masking level and the weapon 

enabling level. 

The output of this step is a table of the available echo excess 

(B..; j= -J, •••» J; k = 1 t K V e 
J K 
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2. 



Scanned Volume Radials; General Case 
In each of the planes j there is now defined for each value of 
range (100k) an available echo excess B^ k for acquisition on the 
azimuth axis. This available excess may be used to find limits for 
off-axis acquisition. Define an angled to be the angle through 
which the weapon turns while a ping travels to a range of 100k and 



returns , 
Then 



*k = 



200 ka)g 



where 

0 ) = weapon search turn rate (deg/sec) 

c = velocity of sound in the medium (ft/sec) 

For echo ranging In the plane j to a range o£ 100k» an energy pulse 

transmitted at angle p will return at the angle p -* k - 

The echo excess A jfc ( p) required to acquire the target at range 100k 

in plane j at angle p off the aximuth axis is 

A k (p) = - 10 log b 3 (p) - 10 logb^CP k ) 
where b 3 ( p) = normalized aximuth pattern function for transmitting 

pattern 

( p) = normalized aximuth pattern function for receiving 
pattern 

If the required echo excess A Jfe ( f» exceeds the available echo excess 
B k then acquisition will not take place at range 100k at an azimuth 
removed p degrees from the reference axis of plane j. However if 
A ( p) does not exceed B jk , acquisition takes place. 
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jn each plane j and at each target range 100k there are two limit ing 
angles between which acquisition will occur with probability one and 
outside of which there will be no acquisition. Looking from the trans- 
ducer toward the potential target define the set of limiting angles 

8 with m 1 for the limiting angle to port and m = 2 for the 
K jmk 

limiting angle to starboard. As shown in Figure 2 f the reference 
axis for polar coordinates in the j th plane is the projection of 
the transducer axis at the time of ping transmission in that plane. 
Positive angles are measured clockwise from the reference axis in all 
planes j and in any horizontal plane. The pattern may not be centered 
on the axis and thusp^ may be positive or negative. Design or 
operational considerations will make it necessary to prohibit reception 
of signals from off the transducer axis in the azimuth plane in excess 
of some limiting angle a in the direction opposite to the search turn 
and « ' in the direction of the search turn. This restricts the 
angles P as follows 

Pjlk “ ^k 

Pj2k “ ^ k | “ * 

Recall that in plane j the maximum range for a non- negative Bj k 
was computed to be 100K^ yards. For a straight running weapon this 
would be the maximum acquisition range in plane j. Actually an 
acquisition range of 100K' yards will not be realized in plane j in 
this problem because the weapon is turning at some rate. The energy 
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100k 



O o 
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Figure 2: Acquisition limiting angles in the plane j „ Line OP is the reference axis 

The x axis is the transducer axis projection,, In this case j2k is zero 



from a given source direction will never return as an echo in the same 
direction. Thus there will always be some penalty for off~axis trans- 
mission or reception. 

Unless the turn rate is extremely high> the insonified volume of 
adjacent weapon pings will overlap. If a given ping scans part of the 
volume insonified by the previous ping> the overlap cannot be credited 
to the volume scanned by the later ping because the overlap volume was 
removed from the volume of target uncertainty by the previous ping* 

The concept of increased probability of detection because of overlap 
is not valid in this case because of the use of the definite range 
law. Thus all overlap volume should be credited to the ping which 
first scanned it. This would give all pings subsequent to the first 
ping during the first search turn an equal scanned volume less than 
that of the first ping. In this paper the situation is approximated 
by dividing the volume of overlap between the pings concerned to give 
all first turn pings an equal volume of insonif ication , This approxima° 
tion decreases P (l) but the error in P (i) decreases with increasing 

n n 

I until the end of the first search turn when it becomes negligible. 

The plane of maximum overlap is plane 0. In this plane^ for small 
k, f*Q 2 k for P in S i will be larger than b + p for Pi n S (i+l). 

See Figure 3* Here b is the angle between the axes of pings i and 
(i+l) or the angle through which the weapon turns between pings. As k 
becomes larger the differences between P ^ for ping i andb + p^^ f° r 
ping (i+l) become smaller. Define p f to be the angle at which the 
two are equal. Thus the right-most angle for acquisition in any plane 
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Figure 3: Overlap of insonified volume between adjacent pings. 



35 



is not greater than p' and the leftmost angle is not less than ?'-*>• 
Some of the extreme planes j may not be useful for acquisition at all. 
By including for consideration only those planes in which acquisition 
may take place the index j takes on the new range from J 1 below the 
axis to J 1 1 above the axis. 

The range or radial distance associated with the angle pj mk is 
100k yards. See Figure 4. There is a range maximum called the ping 
interval range defined by the time between pings (l/PRR,) and the 
velocity of sound in the medium (c). Further, for planes inclined 
above the horizontal there is another range maximum defined by the 
range at which the radial at angle Pjmk breaks the surface of the water 
Similarly for planes inclined below the horizontal there is a range 



maximum defined by the range at which the radial at angle pj mk breaks 
the maximum target depth (D^). Volume scanned below the maximum 
target depth does not contribute to acquisition probability. Thus the 
acquisition range at any angle Pjmk is the minimum of the ranges just 
discussed to be denoted 

Since the index k is no longer of value to indicate range it is 

useful now to relabel the radials such that the index k indicates the 

separation from the leftmost radial in the plane. See Figure 5 . The 

leftmost radial is labeled (R^, Pjl ). The general radial is 

(R B, )• The rightmost radial is (R. r , )• For P lanes close 

' jk 7 r jk J j J j 

to the axis, probably 



PjK. - p ' and p jK _ ~ Pjj. - b 
This may not be true of the extreme planes. 
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Figure 4: Acquisition range at angle Line OR has length 

100k. Line SS 1 denotes the surface. Length 01 is the 
ping interval range. Lined area is scanned. 
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Figure 5: Relabeled radials in plane j„ 
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The output of this step is tables of 
|J k and j = j' > 2 ;k = l, “ * ° > 

plus the angles p and b • 



K. 

J 
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3 . Scanned Volume: General Case 

The radials computed in step 2 define areas of insonif icat ion in 
each of a set of planes j. The area in each plane may be computed as 
the sum of a set of triangles as shown in Figure 6. 

y 1 

Aj = (1/2) pp R jk R j(k+1) Sin (p j(k+l) “ P jk } 



The irregular shape of the areas in the planes j makes it difficult 
to work with them in their present form. Define a set of isosceles 
triangles (one for each plane j) such that the triangle areas are 
equal to the scanned areas A in planes j and the common apex angle 



is the angle turned by the weapon per ping (b). 
altitude of each triangle is 

f "l 1 / 2 



H 



3 



tan(b/2) 



See Figure 7. The 



The set of triangles so defined approximates the areas scanned and 
makes it possible to obtain the total scanned volume SV Q1 on ping zero 
on the first turn. Ping zero is a ficticious ping denoting the 
general case. The computation of SV^ (scanned volume on i th ping of 
first search turn) is identical to that for SV Ql to be developed next. 

Define: p = weapon angle of attack (pitch) (deg). 

The plane j makes an angle of (j+p) degrees with the horizontal. 

The set of isosceles triangles is used to define a set of pyramids. 
A vertical cross-section of the set of pyramids is shown in Figure 8 
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Figure 6: Area of insonif ication in plane j „ 



hi 




Figure 7: 



Equivalence of area and isosceles triangle 
OHH 1 of altitude and apex angle b . 
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7 




ping zero, x is the horizontal axis, y is the 



vertical axis, and OQ is the transducer axis 
extension. 



^3 



and a single pyramid is shown in Figure 9* 
pyramid is 

V. = 1/3 l”area of base altitude 



The volume of a single 



As shown in Figure 9> 2X is the average value of one of the 
base dimensions and Y is the value of the other dimension. 

Thus the base area is 2XY. 



X = |~(Hj + tan(b/2)]/2 

Since sintf> = h/faj = sin 1 )/l 



the altitude is 

h = (Hj sin 1° )A 

and 

Vj = (l/3)(2XY)h = (l/3)(H j+ H (j+1) ) H (j+l) sin 1® tan(&/2) 



The volume scanned is 



J"-l 



SV rt , = 



’01 



j=J' 



The output of this step is tables of H^j j = J * 9 •••, J* 1 
and the volume SVq^. 



kk 
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Figure 9s Pyramid volume V . c Line OE is 0 Line OH 11 is H ( 



4. Special Volume Elements: General Case 

For the case of search at constant depth in a continuous circle 
as is being considered in this paper, the situation changes drastically 
upon completion of the first search turn circle 0 Subsequent search 
turns will effectively scan much less volume because the weapon retraces 
its path of the previous turn. Define the maximum available number of 
search turns (U) to be 



U = 




W OSD 



) PRR /M 
s 




+ 1 



where M is the number of ping cycles in one search turn 



M = ||36o/&|| +1 



IH' 



= greatest integer _ x 



x|| == greatest integer < x 

1 

W = maximum allowable weapon fuel (lb) 

w^ t = fuel expenditure rate at search depth and speed (lb/sec) 
Wqsd = fuel used in dive to search depth (lb) 



(i/3) Dyg v Qt 



v ( 



OSD sln^g 



l/3 D uc = weapon search depth (yd) 

W b 

= average weapon fuel expenditure rate in dive to search 






<P 



DS 



depth (lb/sec) 

= weapon speed in dive to search depth (yd/sec) 
= angle of dive to search depth (deg) 
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Consider a side view or vertical cross-section of the scanned volume 
of a single ping as in Figure 10c The vertically lined area is that 
scanned by the weapon on the i th ping of the first search turn* 

Recall assumption 6 that the target opens radially from the weapon 
position and assumption 12 about weapon orientations on different 
search turns. Define a distance of TES yards as the distance traveled 
by the target during the time it takes the weapon to complete one 
search turn. The weapon scanned area on the i th ping of the second 
search turn is equivalent to that of the areas lined sloping downward 
to the right on Figure 10. The weapon actually scans the same 
geographic volume on each search turn but since the target is assumed 
to open radially the situation is equivalent to translating the 
area scanned on the first turn a distance TES to the left for each 
subsequent turn. If the target is in the area lined sloping downward 
to the left on the i th ping of the first search turn p he will escape 
prior to the i th ping of the second search turn. Similarly for sub- 
sequent search turns. A target above the critical depth CDA or below 
the critical depth CDB may escape but a target at any intermediate depth 
will be acquired eventually if fuel supply permits because of the assump- 
tion of an opening target. The search turns after the first are desirable 
because of the possibility of acquiring a target 1 trapped 0 inside the 
insonified volume of the first search turn. Treatment of scanned 
volumes on the second and subsequent turns will differ depending on 
whether the volume is above or below the depths CDA and CDB. Although 
the areas scanned above depth CDA on search turns 1 and 2 are equal in 
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Figure 10; Vertical cross- section of a single ping. Lines xx“ have length TES 
Lines yy c and lines xx # have length L . 



Figure 10, the volumes represented by these areas are not equal* A top 
view of ping i is a pie shaped slice which makes the volume associated 
with the first search turn area larger than that of the corresponding 
second search turn area. Thus each volume element must be calculated 
independently and the volume scanned by the weapon on ping i of search 
turn u (SV^) is the sum of several independently calculated volumes* 

In this step the volume elements in the general case (i-0) are computed* 
These are modified as necessary for each ping orientation in a later 
step. Although Figures 11 and 12 and the development to follow involve 
the volumes above the search depth, the development for volumes below the 
search depth is identical. 

Consider the vertical cross-section of Figure 10* The cross- 

section has been split into several regions or volume layers identified 

on the left margin by their bottom boundary labeled such that the 

bottom boundary of a level (or layer) at or above depth CDA has a 

length of L, inside the first turn scanned region* The planes j are 
a 

considered using the plane (j = TJ) of maximum H. |sin ( j + p) | as a 
reference plane. In Figure 10 below search depth 
H XJ | sin (TJ + p)| = H (TJ+1) |sin (TJ+l+p) | . 



In this case where more than one of the planes j has the same maximum 

| sin (j + p)|, the reference plane is chosen to be that with the 

least Hj cos (j + p). Above search depth, each plane j such that 

H. sin (i + p)<D„„ - CDA and i>-TL defines a level L , . Also each 

j J WS J d 

plane j such that sin ( j+ [ 1 )^ 0 ^= CDA defines a level L^„ 
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Referring to Figure 11, the volume elements LADV^ , MADV^, and 
RADV^ (d - 1, d; u - 2, eee , U) are significant. The subscript d 

indicates the level and the subscript u indicates the search turn 
number. The level L- is the final level above the critical depth CDA 
and = TES . The corresponding volumes below the search depth are 



LBDV du' mDV du> and RBDV du ( d = 1 > R'> u = 2, 



du 



du 



U ). The scanned 



volume in level above the search depth on the u th search turn is 

ASV = LADV , + MADV, + RADV, 

du du du du . 



One or more of these volume increments may be zero as is LADV^ 

in Figure 11. The corresponding scanned volume below the search depth 

is BSV, . 
du 

Referring to Figure 12, 

the volumes AMV^; d = D+l; •••, D* j u = 1, U are significant. 

The corresponding volumes below the search depth are denoted 
BMV^ (d = D + 1, D' ; u = 1, U ) 

For levels L, below depth CDA and above search depth and for the u th 
d 

search turn, AMV^ is the trapped volume of target uncertainty remain- 
ing after the ping cycle is complete. The volume scanned in level 
above the search depth on the u th search turn is thus 



ASV = AMV ,/ . N - AMV , 

du d^u-1; du 

The corresponding scanned volume below the search depth is BSV^. 

In this step all calculations are done for ping zero and the 
quantities computed for ping zero are saved to be modified later as 
appropriate for the i th ping. The total scanned volume for ping 
zero of the u th search turn is 
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(T-P) 
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(T-P) 
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D' D® 

ST ou= 2 “ v du ♦ jr Kv du 

d=l d=l 

The output of this step is the number of pings (m) in a search 
turn, the maximum available number of search turns U, several tables 
of parameters described in appendix Ck used in computing the special 
volumes, and tables of the following special volumes for ping zeros 



LADV , . MADV, , RADV, : d = 1, 
du 7 du 7 du 7 7 




u = 2, 


® • • j u 


LBDV, , MBDV, , RBDV, : d = 1, 
du 7 du 7 du 7 7 


* * ’} 


u « 2 , 


• • ° u 


AMl^du’ d = D + 1, • • *, D 1 ; u = 


l f 


U 




BMV du 7 d * £ + 1, * * *, D* j u = 


1 • • O 

L 9 9 


U 




ASV du ; d - 1, •••, D' j u = 2, 


• • 9 9 u 






BSV du ; d - 1, •••, D' ; u - 2, 


• • • j u 






SV^; u - 2, U 
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5. First turn scanned volume: i th ping 

Throughout the search phase, the weapon is displaced from the 
center of the target volume cylinder a horizontal distance R^ . The base 
radius of the target cylinder is some distance at the time the 

weapon splashes into the water. The radius then increases at a rate 
( S tc) ec * ua ^ to t * ie target cruise speed to a distance R 7 at the time 
of the first weapon ping. Thereafter the radius increases at a rate 

(S c ) equal to the alerted target escape speed and has a value R. 

1 b i 

upon commencement of the i th ping cycle. 

In the situation of Figure 13, a top view of the target cylinder 
of increasing volume with the first few ping orientations and scanned 
volumes in the first search turn, the insonified volume of ping 1 
includes a large volume outside the target cylinder. Such volume does 
not contribute to acquisition probability because there is an assumed 
zero probability that the target is outside. Therefore the total 
insonified volume must be reduced to the useful scanned volume inside 
the target cylinder. In the situation depicted, the amount of dis~ 
carded volume is reduced with each ping until ping 9 when all of the 
insonified volume is within the target volume cylinder. 

The projection of the radial (r^* p^) pl&ne j to the hori=> 

zontal plane through the transducer is the line OP of Figure Ik 
denoted 

where 

R, jk - R jk c ° s > 
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Figure 13: Top view of target cylinder of increasing 

volume with the insonified volumes of several 
sample pings of the first search turn depicted « 
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Figure 14 : The projection of plane j in the horizontal at 

search depth. Line 00' is the plane j reference 
axis projection. Line IP is Line OP is 



1 .1 . 



Line IP' is R. . Line OP* is R* 



i jk ‘ 
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Use of the law of cosines enables determination of the distance 



I P denoted CR^* This distance is compared to the radius (R^) of the 

target cylinder on the i th ping. If the radius (R^) smaller^, it 

defines the new triangle IOP 1 and the projection of the side 0P ff back to the 

plane j defines a new acquisition range (R ,) to be associated with the 

l jk 

angle p., . That is 

J t 

R ijk 



R. 



ijk cos 



ijTTi 



However if the distance CR t1 is the lesser of the two then 

jk 



R. = R., 
ijk jk 



The remaining portion of the computation for scanned volume (SV^^) 
on the i th ping of the first search turn is identical to that in 
step 3 for 

The output of this step is the table of volumes SV^j 1 = 1^ M=X 

and the distances R and R.: i = 1, °°°f M-l. 

o 1 
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6. Scanned volumes subsequent to first turn* i th ping 

Just as the first turn scanned volume had to be corrected for the 1 th 
ping because of the effect of the target cylinder radius^ so must the 
scanned volumes on subsequent turns* 

Define the intersection of the ping i bisector with the target 
cylinder radius at transmission of ping i to be a distance CH from the 
weapon . 

Referring to the vertical cross-section of the insonified volume 
in Figure 15, verticals through four of the possible intersections have 
been plotted and labeled CHI through CH4 . The effect of each of these 
four cases will be discussed as it applies to a level above and a level 
below the critical depth CDA (i,ej levels 2 and 4)* The treatment of 
the situation below the search depth is identical* 

In the case of CHI the first turn scanned volume is reduced* 

However the scanned volumes for all subsequent turns are the same as 
those computed in step four for the general ping* This is because 
CHI lies to the right of all areas representing equivalent volumes* 

In the case of CH2 there is still no reduction for volumes in 
level 2* However there is a small reduction in the volume scanned on 
the second search turn in level 4 , The volume reduction is that 
volume associated with the area in level 4 to the right of CH3o 

It is computed as a reduction in AMV^j (the volume missed in 
level above search depth on the first search turn)* 
ftuiS (ASV 22 ) plng t equals (ASV 22 ) plng 0 but 

(AS V>ping i is Uss than (“Wng O' 
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as numbered o Escape volumes have primed numbers according to the search turn of escape 0 



In the case of CH3 there is a reduction in the scanned volume on 

the second search turn in both levels and 0 In level the volume 

(ASV^) . . is less than (ASV^) . A by an amount equal to 

v 22 'ping i v 22 'ping 0 J M 

RADV 22 plus a portion of MADV^* 

In level the reduction still amounts to a portion of AMV^ 0 

In the case of CRk there is no scanned volume in level above 
search depth on the second turn, That is ASV ’ s 0 for ping i„ 

Further there is a small reduction in level in the scanned volume 
for the third search turn. For level there remains some scanned 
volume on the second turn and there is also a reduction in scanned 
volume for the third turn. 

The output of this step is a table of scanned volumes 
SV. u ; i - 1, •••, M-lj u « 2, U. 
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7- Target Volume: i th ping 

Consider that portion of the target cylinder volume computed at 
a base radius of within the radials bounding the ith ping* See 
Figure 16* The area contained between the boundary radials of the 
i th ping and within the cylinder of radius multiplied by the 
height of the cylinder is this volume * The probability that the 
target is in this i th volume slice of the target cylinder is the 
volume of the slice divided by the total cylinder volume 0 This 
probability (PQ^; i = 1, •••, M-l) remains constant throughout the 
problem because of the opening target assumption number six* 

From Figure 16: 

pq _ area of horizontally lined region 

1 7 T Rj 2 

The target volume (TV^.; i = 1, M-l) associated with the i th 

ping on the first search turn lies between the same radial limits but 
within a target cylinder of radius instead of R^. That is 

TV^ = (area of vertically lined region) 

Figure 16 shows the angle ( b' ) assigned to the M th ping to be less 
thanb . It is 

b ' = 360 - (M-l)b 

This smaller angle is required to complete the search turn so 
that the i th ping of the u th turn is oriented in the same direction 
as the i th ping in the (u + 1) th turn. Thus the M th ping is in 
part ficticious. 

The output of this step is thus the tables of TV^j i = 1, •••, M-l 
and of PQ^; i = 1, •••, M-l and the angle b* . 
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Figure 16; Top view of target cylinder of increasing volume . 

Lines CT and Cl' are boundary radials of ping i. 

2 

Horizontally lined area is PQ^ times TTR^ „ 
Vertically lined area times cylinder height is TV^» 
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8. Acquisition Probabilities: I pings 

The instantaneous probability of acquisition on the i th ping is 
simply the ratio of the volume effectively scanned to the volume of 
target uncertainty. 

Define: g^ = P ^acquisition during the i th ping cycle of the u th 

search turn^j 

Then g lu = SV lu ^ J i = 1, •••, M-i; u = 1, •••, U 

Recall that PQ^ is the probability that the target is within the 
radial limits of ping i. 

The volume TV^ is the total volume of target uncertainty within 
the radial limits of ping i, The volume SV^ is the volume scanned 
on the first turn. The volumes scanned on subsequent turns are volumes 
measured inside the first turn scanned volume and under the conditions 
existing for the first turn. These volumes are then moved intact with 
the opening taget . Thus the applicability to all turns of the volume 
TV^ measured using the first turn cylinder radius. 

The instantaneous probability of detection on the M th ping is 
approximated by 

«H> = *(M-l)u {b,/b) * g l(u»l) (& * b ' )/b 

The M th ping is alloted the same time interval as any other ping 0 
Errors involved should be small. 
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There remains the requirement to relabel the to be 

g^; i = 1, •••, I, to compute I as a function of as described in 
the beginning of the chapter, to compute P as a function of and 

A 

to tabulate the values of versus as is allowed to vary from 

i 

W to W . This final table is the input to the dynamic programming 
model . 
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9. Input Requirements 

A. Medium 

c = velocity of sound in water (yd/sec ) 
a = attenuation coefficient (db/yd) 

10 log n> s = surface reverberation scattering coefficient (db) 
10 log my = volume reverberation scattering coefficient (db) 

B. Target 

TS = target strength (db) 

TN = radiated noise level of the target (db at 1 yd) 

D_. = maximum target depth (ft) 

S = target speed during weapon search (yd/sec ) 

i. j 

S = target cruise speed (prior to becoming alerted on 
the first ping cycle) (yd/sec) 

C . Weapon 

WS = transmitted on-axis source level of weapon (db at 1 yd) 
AT = ping length (sec) 

PRR = weapon active sonar pulse repetition rate tor 
s 

search (pings/sec) 

A f = receiving bandwidth (cycles/sec J 

WEL^ = weapon enabling level tor a range of look yds (dbj 
SN = equivalent isotropic self°noise level at search 
speed, depth, and turn rate (dbj 
RD = recognition differential for 50$ probability of 
acquisition (db) 
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b ^ = normalized pattern function in the vertical plane 

for transmitting pattern (a numeric) 

= normalized pattern function in the vertical plane for 
receiving pattern (a numeric) 

b^ =s normalized pattern function in the azimuth plane for 
transmitting pattern (a numeric) 
b^ = normalized pattern function in the azimuth plane for 
receiving pattern (a numeric) 

J = limiting angle for off=axis reception in the vertical 
plane (deg) 

cr» - limiting angles for off=axis reception in the azimuth 
plane (deg) 

DR = directivity index for receiving (db) 
p = angle of attack (pitch) of weapon (deg) 

0) c = search rate of turn (deg/sec) 

^ DS = angle of dive in diving to search depth (deg) 

D = weapon search depth (ft) 

W j 

= weapon speed in dive to search depth (yd/sec) 

S T ,„ = weapon speed during search 

WS 



W = maximum allowable weapon fuel weight (lb) 
w = weapon fuel expenditure rate at search depth and 
speed (lb/sec) 

w = average weapon fuel expenditure rate in dive to 
search depth (lb/sec) 
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D. Delivery weapons system capability 

Rj = radius of target volume cylinder at time of weapon splash 
Statistical averages and distributions are available for the medium 
inputs. Intelligence reports and estimates or CNO requirements may 
provide the target inputs except for the target speeds which may be 
selected in the higher ranges in accordance with criterion one 0 
The weapon inputs are available from weapon designers at least in the 
form of educated estimates. Finally the input may be estimated 
from observations of the performance of available delivery weapons 
systems . 
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CHAPTER IV 



PROBABILITY OF HIT 

Consider a tactical situation in which a weapon acquires a submarine 
target at some initial slant range and make the following assumptions * 
Assumption 1: The weapon makes an attack and a series of independent 
reattacks on the target culminating in either a hit or fuel exhaustion* 
Assumption 2: At the time of acquisition, the target is opening at 

maximum speed. Target course and speed are constant throughout the 
attack and reattack. 

Assumption 3: On the attack and after reacquisition on each reattack 

the weapon maintains contact and closes at attack (maximum) speed* 
Assumption k: The target operates in depth within specified minimum 

and maximum limits such that the maximum depth limit is not greater 
than the weapon maximum operating depth* 

Assumption 5: The depth at which a hit may take place is a random 

variable denoted the hit depth which has an assumed probability 
density function. 

Assumption 6: The average weapon depth for the attack is the arithmetic 

average of the search depth and the average hit depth* The weapon and 
target depths for the final stage of the attack and throughout each 
reattack are equal at some depth within the limits of assumption four* 
Assumption 7 : The attack and each reattack is from astern the target* 

There is some close-in range at which the weapon sensors lose contact 
with the target. On losing contact with the target the weapon continues 
on a straight track at attack speed for a preset dead time* Upon com- 
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pletion of the dead time run, the weapon instantaneously slows to 
search speed and commences a left turn at reattack turn rate and 
constant depth. 

Assumption8 : The weapon reattack turn is a circle with no advance or 

transfer. The weapon reacquires the target upon completion of one 
complete reattack turn and instantaneously increases to attack speed. 
Assumption 9: At the time the close-in weapon sensors lose contact 

with the target, there exists an intersection of the weapon track 
extension with the vertical plane through the submarine target center 
which is perpendicular to any plane containing both the target and the 
weapon. Let the target projection in the vertical plane be approximated 
by a circle denoted the target circle centered on the target longitudinal 
axis with a diameter equal to the length of the beam dimension of the 
target. Let the circle center be the origin of a doubly infinite 
axis containing the intersection. Assume that the intersection co- 
ordinate on the axis is a normally distributed random variable with 
mean zero and standard deviation o «> 

Assumption 10: A hit occurs if the weapon track at the closest point 

of approach to the target center passes within a distance of the center 
not greater than the radius of the target circle plus a small proximity 
correction factor. 

Define : 
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The derivation is analogous to that for acquisition probability 0 
Define: 



t = time the weapon sensors lose contact with the target on an 
attack or reattack (zero reference time) 

RH = range from the weapon to the target center at time t Q (yd) 
t = time from t of closest point of approach of weapon and 

L»iA O 

target (sec) 

MRH = distance traveled by weapon from time t Q to time t^p ( yd ) 

S = weapon attack speed (yd/sec) 

WA 

S_, = target maximum speed (yd/sec) 

BW = target beam width (yd) 

VT = weapon proximity correction factor (yd) 

r = coordinate of intersection of the weapon track extension 

with the vertical plane (A) through the submarine target center 
perpendicular to any plane containing both target and weapon at 
time t Q (yd) 

r* = coordinate of intersection of the weapon with the vertical plane 
(A* ) parallel to plane A containing both the submarine target 
center and weapon at time t^p (yd) 

F^(x) = cumulative normal probability distribution with mean zero and 
standard deviation one evaluated at the point x 
Consider the attack situation depicted in Figure 17« 

The xy plane is the vertical containing the weapon and target at time 
t Q * Planes A and A f are parallel to the yz plane and contain the 
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Figure 17 % Planes A and A 8 in the attack situation 



bases of two right circular cones with a common apex at the origin 0 The 
line OR is the weapon track and its orientation with respect to the x 
axis determines the orientation of the axes LL and L 8 L 0 on which 
coordinates r and r 1 are measured. The common intersection of the cone 
of base radius r* 9 the plane A 1 9 and the xy plane is a line of length 

2 | r 1 1 . By similar triangles 

i MRH 
r = RH r 

Since 



MRH = RH S„ A / (S„ A - S TM ) 



then 



r ' = r S WA / ^ S WA " S TM^ 



Thus r 1 is a random variable with mean zero and standard deviation 



“ 0 S WA / ^ S WA “ S TM^ 



Therefore 



h = 



P 2 | r* | < BW + 2 Vt] 

= P |r' | £ (BW/2) + VT] 



= 2 F 



0 



(BW/2) + VT 
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h 



n 



2 F 






(BW/2) + VT 
0 S WA / ^ S WA " S TM^ 



- 1 



Define: 

HD = hit depth (ft) 

^HD^ X ^ = probability density function of the random 
variable hit depth 
D^ = maximum target depth (ft) 

D MIN = minimum target depth to be considered (ft) 

Define a set of Q equal depth subintervals in the depth interval 
D MIN* °TmJ suc ^ fc hat the depth at the lower boundary of the q th 

°MIN + q D TM~ D MINj ‘ 

= P hit depth is in the q th depth subintervafj j q-1^ ° ° ° > Q 
Assume as in Figure 18: 



subinterval is 



Define 



~ 2 (x - Dwnj) / (D- 



fc HD w - - - -MIN 7 ' VU TM " D im y 



Then; 




S3 - l . 
Q 2 ’ 



q 



i. 



Q 



Define : 

w' = fuel used in the n th attack or reattack (lb) 
2n ' 

=s initial acquisition range (yd) 
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f HD^ x ' 

A 



2 

d tm" d min 



0 



Figure 18 : 




Assumed probability density function of the 
random variable hit depth. 
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w^ t = weapon fuel expenditure rate at attack speed and 
at the average depth for the attack, (ib/sec) 

HD = average hit depth (ft) 

Then 

5B = 2/3 D ra ♦ 1/3 d min 

and the average depth for the attack is 



WS 



+ HD 



- 1/2 D wg + 1/3 D tm + 1/6 D. 



TM 



MIN 



For the attack 



w 



W = 
21 



2t 



WA 



TM 



Define: 



w 



Stq 



w 



Atq 



'DT 



'CT 






= weapon fuel expenditure rate at search speed and at 
the midpoint of the q th depth interval . (lb/sec) 

= weapon fuel expenditure rate at attack speed and 
at the midpoint of the q th depth interval, (lb/sec) 

= dead time from t at which the weapon commences the 
reattack turn (sec) 

= time from t at which the weapon completes the reattack 
turn and reacquires the target (sec) 

= reattack rate of turn (deg/sec) 



Then; 



t 



CT 



'DT 



+ 360° 



t 



CPA 



RH 

S WA ' S TM 
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The weapon will spend the following time at attack speed during each 
reattack: See Figure 19 . 

^ T CT " fc CPA^ S TM~ ^ t DT° t CPA \ S WA 
A = DT - CPA + S WA - 

The weapon will spend the following time at search speed during 
each reattack: 

T - 360° 

S ^ 

Therefore : 



Q 

W 2n * 21 [ T A W Atq + T S “stq ] P Eq 
q=l 



n = 2, •••, N 
Define : 

S = minimum preset value of attack speed (yd/sec) 
o 

w^ = available fuel weight for hit, event 2, (lb) 

w” = available weight for modifications to the propulsion 

system to provide a maximum speed capability in excess 

of S (lb) 
o 

N is a function of w^ according to 



N 



N+l 



w_ 



2n 



— << 



w, 



2n 



n=0 



n=0 



However, recall that w^ n is a function of attack speed which is it= 
self a function of w^' . 
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Figure 19: Geographic plot of weapon and target tracks 

separated to show time differences • The 
distancesbelow the target track are measured 
from the joint position at time t . 

v iA 
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Thus 



P r (N) = 1 



N(v£, w£' ) 



n=l 



t 



h n (wj') 



Recall w g = + w” 

I _l 

W> W = lower and upper limits respectively on the range of fuel 
weight (lb) 

n Jl 

W , W = lower and upper limits respectively on the range of 
weight for modifications to the propulsion system to 
provide a maximum speed capability in excess of S q (lb) 
For each value in the range for w^ 
min -£w' , W' <w 2 <~min W + W , 

the function P is evaluated for each possible combination 
H 

of the variables w^ and w^ f as they are allowed to vary such that 

they remain within their respective ranges 
1 1 

W ^ w_ ^ W 
- — 2 

it ii t 

W — w 2 

and their sum is equal to w^* Let K be the number of values of 
w^. This procedure provides sets of values 



w 



2jk* W 2jk* *Hjk 
for the k th value of w^* (k = 1, •••, K) 
The entire set of sets of values 



> j “ ^y * * ° y j 



W 2k* W 2jk" W 2jk> P Hjk j j = 1# 1 J k 

k = 1, 4 • *, K 
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is the input to the dynamic programming allocation model 0 

Input Requirements 

1. Target 

BW = target beam width (yd) 

D TM = tar g e t maximum depth (ft) 

= target maximum speed (yd/sec) 

2 . Weapon 

S = weapon search speed (yd/sec) 
w j 

VT = weapon proximity correction factor (yd) 

00^ = weapon reattack turn rate (deg/sec) 

i _i 

W, W = lower and upper limits respectively on the allowable 
range of fuel weight (lb) 

S q = minimum allowable weapon attack speed (yd/sec) 

ii _i i 

W , W = lower and upper limits respectively on the range 
of weight for modifications to the propulsion 
system to provide a maximum speed capability in 
excess of S (lb) 

G^ = Relationship between increased weapon speed above S q and 
increased weapon weight 

= Relationship between fuel expenditure rate and speed and 
depth. Recall that this relationship may vary as a 
function of attack speed only in the range for speeds 
above search speed. This relationship is required to 
specify and w 2j . (lb/sec) 
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3. Miscellaneous 



RH = range from the weapon to the target center at the instant 



contact with the target (yd) 

0 = standard deviation of the normally distributed random 

variable: coordinate of intersection of the weapon track 
extension with the vertical plane (A) through the 
submarine target center perpendicular to any plane con- 
taining both target and weapon at time t • (yd) 

D llTKT = minimum target depth to be considered (ft) 

MIN 

= initial acquisition range (yd) 

Q = number of depth subintervals in the depth interval 



Intelligence reports and estimates or CNO requirements may 
provide the target inputs. The weapon inputs are available from 
weapon designers at least in the form of educated estimates. The 
inputs RG and O are functions of the medium, the weapon sonar and 
guidance systems, the target physical dimensions, and the orientation 
relationships between the target and weapon. In the absence of a 
specific theory these must be handled by simulation or sensitivity 
analysis. The input is a function of the capability of the 

weapon delivery system and the tactical situation as seen by the 
target commander. In the absence of specific considerations cri- 
terion one would determine to be zero. The input R A should be 

MIN A 



(t Q ) during close-in attack when the weapon sensors lose 
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selected close to the maximum weapon capability in accordance with 

criterion one* The function P may be very sensitive to changes in 

R and an examination of results over a range of values for R may be 
A A 

appropriate* The input Q should be large enough to point up the 
effect of an irregularly shaped fuel expenditure rate curve but small 
enough not to excessively prolong the computer program* 
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CHAPTER V 



PROBABILITY OF KILL DAMAGE 

Consider a tactical situation in which a weapon explodes within 
a specified proximity of some portion of a submarine target and make 
the following assumptions* 

Assumption 1: The hit depth is a random variable which has an 

assumed probability density function* 

Assumption 2: The point on the target at which the hit occurs de- 

noted the hit point is a random variable which has a uniform distri- 
bution over the entire surface area of the target* The weapon 
approaches the target from astern, loses contact at some close-in 
range and continues on a straight track* Since the target may be 
in a turn, climb or dive, all points of the target are available 
as hit points * 

Assumption 3: The target operates in depth within specified minimum 

and maximum limits such that the maximum depth limit is not greater 
than the weapon maximum operating depth* Divide the target surface 
area into P area elements, each of which has a constant hull 
vulnerability throughout. Divide the interval of possible hit 

depths into Q equal depth subintervals such that the depth at the 

i 

lower boundary of the q th subinterval is 



d min + q 



d tm" d min 



82 



where 



D TM = maximum target depth (ft) 



MIN 



= minimum target depth to be considered (ft) 



Define: (w^) = P £kill damage given an explosive weight wjj 

f„ (x) = probability density function of the random 
HD 

variable hit depth 

P £ ^ = P £hit depth is in the q th depth subinterval J | 

q = i> Q 

P^ = P j^hit point is in the p th surface area elementj; 



P = 1> 



P 



= hull vulnerability for kill damage of the p th 
target surface area element at the midpoint depth 
of the q th depth subinterval given an explosive 
weight w^. This number is one if kill damage occurs 
under the specified conditions and is zero if kill 
damage does not occur • 
w^ = explosive weight (lb) 

It! V « • 

W , W = lower and upper limits respectively on the range 
of explosive weight (lb) 

Assuming as in Chapter IV 



£ hd ( x) - 



2 ( *- D mt l ) 
(d tm* d min> ! 
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Then 



P = — *1 — - • a = 1 ••• o 

% q 2 ' q h ’ * 

p _ area of p th surface area element ; p s 1, P 

Ep “ total target surface area 



P E = P E ^ w 3 ^ = 



HV (w.) P_ P_ 
pq 3 Ep Eq 



p=i q=i 



I M M 



The table of P versus w as w varies from W to W 
b 5 5 — 

is the input to the dynamic programming model* 



Input Requirements 

1. Target 

P E p = P j^hit point is in the p th surface area element] 
p = 1, • • •, P 



= maximum target depth (ft) 

TM 

2. Weapon 

iii _i i • 

W , W = lower and upper limits respectively on the range 
of explosive weight (lb) 

3o Target and Weapon 

HV (w^) = hull vulnerability for kill damage of the p th 

pq 3 ' 

target area element at the midpoint of the q th 
depth subinterval given an explosive weight (w^) 

4* Miscellaneous 

D = minimum target depth to be considered (ft) 

MIN 

Intelligence or CNO requirements may provide D_ w . P and some 

TM Ep 

of the information for HV (w^) . Weapon designers may estimate 

it* _i i i pq j 

W and W and weapons engineers may provide extrapolated 
estimates of HV (w_). 

pq 3 
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APPENDIX A 



FLOW CHART NOTATION 



The following notation is used in the flow charts of the remaining 
appendices . 



Function Box 



A — 



The function box contains a number or expression to be evaluated 
A which is to be placed in memory cell Bo 

Choice Box 




The choice box contains a comparison of the number or express 
sion A to the number or expression B, Exits are chosen as follows s 

Exit is chosen if 
A> B 
A > B 
A = B 
A <B 
A < B 



Exit Label 



< 

< 
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In Box 



INPUT 



The in box contains an input to the program * In appendices 
C 1 through C 8 the output from one appendix is carried over intact 
into the next appendix. Inputs to be used for the first time in the 
sequence of flow charts are placed in a separate in box* 

Out Box 

OUTPUT 



The out box, one for each appendix, contains the output of the 
appendix program. 



Index Box 



{ 



1 + i 



} 



The index box contains a manipulation on an integer in an 
index register and is the tool which provides for loops in the 
program . 
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Subroutine Box 




The subroutine box indicates a return jump to the subroutine as 
numbered by k. The flow chart of each subroutine is shown after 
the flow chart of the first routine using it , 

Return Box 




A return box indicates the end of a subroutine implying a jump 
back to the main routine. 

Fixed Connector 





A fixed connector is used to connect remote parts of the flow 
chart with one another. 



Assertion Box 



Note j 



1 

An assertion box refers the reader to a note found at the end of 



the flow chart. 
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Instruction flow is generally from left to right and from top to 
bottom. All boxes are entered from the left or from the top and all 
exits from boxes are made down or to the right. The arrow is used to 
further indicate the sequence of computer instructions. Arrows are 
used at all line intersections to indicate entry into the program 
between two boxes. 
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APPENDIX B 

THE ALLOCATION MDDEL 
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Note 1: Provision is made for a set of payload weight limits (W ) in 

n 

case it is desired to examine a range of values. The input 

W, is to be the greatest of the W . 

N n 

Note 2: The subroutine RJl computes the h^(w) and the w„(w) for 

i=l, 3. h^ is handled separately. 
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Note 3 



Note Ij- 

Note 5 
Note 6 

Note 7 
Note 8 

Note 9 



This loop provides that, for the W considered, the function tr 

is maximized over the w^ . For h^ the maximum meaningful value 

in both loops is w^=D so th e subroutine is stopped by the 

comparison W : D. However for h^, the maximum meaningful 

value in the two loops is different, the one in this loop being 

w-._ and the one in the other loop being W...„ s D. 

3 b MAX 

This loop stores the h^(w) and w^(w) obtained as W is allowed 
to vary over its range of meaningful values, 

This loop maximizes the w^^ over j for a given k„ 

The maximum meaningful value of w^ for a given W is 



min 






The h^(w), w^Cw), and w^’Cw) are stored for each value of W. 

This loop works backwards through the stages to retrieve 

the optimum weight values for each input payload maximum 0 

The are the optimum volume sums to be compared with the 

input payload maximum allowed volume V 0 If the solution is 

not optimal with respect to the restraints^ a new X is chosen 

and the program repeated. If the final solution is obtained^ 

then fuel weight is + 1^, weight for added attack speed is 

H and explosive weight is P « 
n n 
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APPENDIX Cl 

AVAILABLE ECHO EXCESS (p^) 
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Note 1; This flow chart neglects the effects of ambient sea noise and 
of volume reverberation. Both are considered to have negli- 
gible effect on the interference background masking level. 

If it is desired to consider the volume reverberation, it is 
M (R v )k = WS - 21^ + 10 log m v “(J v ) k + 10 log 

+ 20 log (100k) + RD; k 1, •«», K 

where M(R ) = volume reverberation masking level at range 100k yards (db) 

' v'k 
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WS = on~axis source level (db at 1 yard from transducer) 

= one way transmission loss at range 100k yards (db) 

10 log m^ = volume scattering coefficient (db) 

(j )^ = volume reverberation index of the transducer (db) 
c = velocity of sound in the medium (yd/sec) 

AT = ping length (sec) 

RD = recognition differential of transducer (db) 



The index (j y )^ may be approximated by 



= 200 ko) 



( J v ) k = " 10 lo s b 3 \ ( P -s }d _n_ 



k = l. 



K 



where m = search rate of turn (deg/sec) 

j = off-axis angle in the vertical plane (deg) 

P = off-axis angle in the azimuth plane at time of 
pulse transmission (deg) 

b^, ^ 2 * s normalized pattern functions of 

the transmitting (l, 3) and receiving 
(2<> 4) patterns in the vertical (1,2) 
and azimuth (3$ 4) planes 

_Q_ = solid angle subtended by the sound energy beam 0 

This in effect assumes a medium of infinite extent 
in the direction of the sound beam 0 
Note 2s M(SN) = SN + 10 log A f-DR + RD 

where M(SN) = self-noise masking level (db) 

SN = equivalent isotropic self-noise level at search speed 
and depth (db) 
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Note 3t 
where 



Note kt 
where 

Note 5: 

Note 6: 
where 

Note 7: 
where 

Note 8: 



£ = receiving bandwidth (cycles/sec) 
DR = receiving directivity index (db) 
2H « WS + TS-M(SN) 



H = maximum possible on^axis range for acquisition (yd). 

This neglects noise levels except the self=noise level 
and neglects the weapon enabling level. 

TS = target strength (db) 

= 20 log 100k + 100 ak 

a - attenuation coefficient (db/yard) and it is assumed that 
the transmission anomaly is zero. 

This loop finds the upper limit (lOOK) on the acquisition 
ranges to be considered. 

M(TN ) k = TN + 10 log A f-H k + RDj k-1, •••, K 
M(TN) k = target noise masking level (db) 

TN = radiated noise level of the target (db at 1 yard from 
transducer) 



<#> . = arc sin 
k 



D WS 

300 k 



= the angle above the horizontal at which an arc of radius 

100k yards and of origin at search depth will inter= 

sect the surface of the water 

D „ = weapon search depth (ft) 

W o 

(j ) = surface reverberation index of the transducer (db) 

\i “angle of attack (pitch) of the weapon (deg) 

This is an approximation described in Chapter k- of refer- 
ence k . 
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Note 9: M(rs)^ = surface reverberation masking level at a range of 

100k yards (db) 

10 log m^ = surface reverberation scattering coefficient (db) 
Note 10: This subroutine adds the effects of two noise levels in 

decibels. A rounding off process is employed based on Figure 3 
of reference 3* A difference in levels of more than 10 db pro- 
vides a sum effect equal to the larger of the individual 
levels. Two equal levels have a sum effect equal to 3 db 
higher than either individual level. 

Note 11; WEL^ = weapon enabling level (db). Values out to a range 

of 100K 1 yards are inputs where is estimated such 
that it will be larger than the K found in the program. 
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APPENDIX C2 

SCANNED VOLUME RADIALSp GENERAL CASE 




INPUT 



PRR fl , « 9 « ' p D, 



INPUT 



TM 



APPENDIX Cl OUTPUT j aytRR^-^b 



200 a^/c { 1 > 1 y ~ = Kq -c- (i => 1) 0^ — 



P 



Note 1 






- ( 1 ♦ i-^i> 



_T> 



i - 1^1 



> ~»J> 



f 



Note 2 



(RJ3) <l->k 




/ 1 » k ->¥V 



P* 1 <1-»1 >-%’ - (llb^ll ♦ 1)^ ♦ (i - l)« l -^p 1 



fT~+ 1->Q - 







J -3> J* — j — J-3-J* 



i 





L 



? Jk 9 jk s 

f>’. b, 9 , <D S , Dv|s , D m , PI®, 
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RJ4) 



1 — in - ' 1 ^ k ) — ^ 0 ja x 1 



i 



p jmk : p Jm(k+l) 



j'Note 9 1 



(~1 ± k ->k ) — ( 1 » n — 5>n / — 



P 



Note 10 



Pj a k-^P J (k-n)h L - | 100k • c/2PRR a p^ c/2PRR s -^R j(k ^ } 




100k — ->■ R 



'J (k-n) 



J(k=*n) • -sin(j+n)ees 




(d tm - D ws } / 3 



R, 



sin ( j +ji)cos^ k ^ D ^ J (k=m) 



V s 

P • 


> 


“vt/ 3 


^JUc-n) sin(J+|i)cosp^ ktta ^ 




sin(J^)cospj ~^ R j (k»n) 




[Note 12 



p j(k-n) • p ’ 




< 






( 1 ♦ k 


:^k) 



k - n — ^K„ 




(RTN 



Pjl(k-n) * p J2k 



< k - 1 -5>k^ ( n 2 — >a) (J) 







- 1 — $»n) (2 -g» m ) - 



Note 13 



[k - l^=>k~) (2 -»a ) (n ° 2 — >nV 



<D 



100 





which meet the requirement that the echo will be received 
in the interval CC *J • 

Note 2: The overlap in plane zero is computed first to provide the 

angle p* . Then p* <= b and p® define the limiting angles 

of the P jmk" 

Note 3: A Jki ■ A Jk ( P ) for the i th P . 
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Note 

Note 

Note 

Note 

Note 

Note 

Note 

Note 



4; After the computation of which indicates that the maximum 

acquisition range in plane j and all of the p in plane j 

jmk 

have been computed^ the subroutine ends 0 
5s This indicates that there is no acquisition probability in the 
plane j and changes the range of planes to be considered from 
-J* ° • • 9 J to j\ • ° ° 9 J J 5 . 

6s In this case the p is limited by the P c and would have 

been larger if the overlap between pings had not been removed 

from the considered scanned area 0 

7: This loop starts with the leftmost angle ( p 0 ®(l«l) ^ ) 

to be considered and computes each angle to be considered 

indexed by i and incremented by 0^ 9 to a maximum angle p T s P° © 

8: This point is the commencement of the loop which will compute 

the p and R starting with j * =J 0 
jk jk 

9: There may be several equal P because the scan azimuth limits 

were reduced by the overlap 0 The first of the equal p ^ to be 
considered will then be the one of maximum k 0 
10; The ping interval range is - yards /ping. 

11; The range at which the radial at angle p ^ intersects the 
surface of the water for j >• = }& is computed below. 

See Figure 20. The distance 0C is desired. The search depth 
is (DWS/3) yards. 



length (CB) = Dyg / 3 sin(j^p) 
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y 




Figure 20. Computation of radial length for angle p 

The plane j for which the scan outline is shown 

is perpendicular to the vertical xy plane. Angles 

CBO and DBO are right angles . Length Sis tefee 

search depth of D yjg/3 yards. 

I 
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Note 



Note 



length (OC) = length (0B)/eosp jk = Dyg/3 sin(jj-frji) eosp 



For j - the development is the same except that (j + p* ) 

is negative and the considered depth is the difference between 

the maximum target depth and the weapon search depth 0 

TM 

12 s If to the left of the radial of maximum range K. then m - 1 

J 

and the following angle will have a larger k c 
If to the right of the radial of maximum range K. then m ^ 2 
and the following angle will have a lesser k 0 There will be 
one or two radials of maximum range 0 That is 



PjlKj — ^J2Kj 



13 : If Pj lK = Pj 2R » the program 



program includes it, 



8 o Otherwise the 

j2k j 
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APPENDIX C3 



SCANNED VOLUMES^ GENERAL CASE 




Note Is SVq^ i- s not required in the program but it may be desired 
to print it here as information. 
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APPENDIX C4 



SPECIAL VOLUME ELEMENTS l GENERAL CASE 
INPUT INPUT 



£H 



> K » S 



WD‘ 



^DS ! 



W» 






Lib 



APPENDIX C3 OUTPUT 



} 



^ 11360/^)11 



3 



+ 1 — ^>M 




* x 5*^ n 





D WS W Ol/ 3 8in ^DS S WD~ >M OSD 



Cr 



[(w , “W nsn )/w 



OSD 7 ' It 



PRR g /M]| ♦ 1-^-U 



-TES 






]D- 



t£ 



» 0 







RJ6 



D« 




OUTPUT 

^du’ ^du’ 
LBDV, 

AM d 

BM7 d 

ASV r 



du v du 9 du v 9 

^du 9 MBDV du9 RBDV du ; d=l/ 

7 du j d = d+i 9 ooo 9 d»; u = i 9 

7 du J d = D+l 9 D'S u = 1, 

O 

T • J ^ 1 006 nil © ^ 




L 



“"W* 


= 1 , 


o o o r\ o @ 

9 U $ 


u : 


s 2 9 


ooo 9 u 


qSV 0u J q = 


= a 9 


B 9 u = 


0 0 


% U 




h Ad 9 


B AId 


9 tan ^AId 9 d = 


1 . 


0 0 © 

9 


D* 


B A0d 


, w AM ; 


d - 1 9 0 


0 0 

9 


D 




h Bd 9 


B BId 


9 tan^g 9 d 


1 , 


000 

9 


D« 


B B0d 


* W B 0 d 5 


d - 1 9 0 


0 0 

9 


D 

c=> 




P Jk 9 


R jk 9 


j = 


J«, 00 % 


J» « 


0 k s 
9 K 


“ij 0 0 © 

x 9 9 


^DS 9 


S WD 9 


S TS 9 


w 0SD 9 w lt 9 


m 9 u 




P*» 


^9 


D ws 9 




s 







j 
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H (j-ATJ-d+ 3 ) 3ln ^“ ATJ “ d+s+ ‘ i ^ 

H (J-ATJ-d +S ) COS 

H (j«d-r) sln ^ 4d ~ r+ ^~^ SS d 

H (J+d-r) C0S ^ +d ” r+ » A ^ sc d 

H ( J>d.r-l) sln SS k 

H(j^ r _i) co s(J+d-r-lnO-5^ SC* 



Not© 7 
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L-^RTN) 
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Note 10 
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Note 13 



RJ13 — L - 



A 2 tan(6/2) h q(J -^B 




(RTN 




RJ14H- 



I'Note 14 



(1/3) tan (6/2) C(D 3 - E 3 )^=>F 




RTN 




Note 11 



(RJ17) — 



(TES =• L ( d= 2)) = BS d ) 

(d-1) “ 



HL 



HL, 



— d-^D 



{YTT^rJ— (HL d - BS d )/(Hl (d _ 1) - BS d ) ^»CT 



~ ' HL (d-l) * HL d ^‘ h qd 




(RTN 




depth. 

Note 2: The following comparison is used to determine the plane 

j to be used as a reference plane for the computations. 



Ill 







Note 3t 



Note k 



Note 5 



Note 6 



Note 7 



The following four function boxes set up reference constants 
for the subsequent computations <> 

: The index d is used to identify the levels 0 The index (d«r) 

identifies the current plane j above the plane TJ« The 
index (d-s) indentifies the current plane j below the plane 
TJ- ATJo 



i This comparison determines if the level d is to be determined 
by a plane j > TJ or by a plane j <(TJ= AtL) 0 

It is only required for depths at or above the critical depths 
that is for d ^ D . 

: This subroutine computes the distance L^ 0 See Figure 21 • 

As shown H^ tj + ^ sin (TJ + 2 + p ) > 

H(tj g) s ^ n + p ) and ATJ - 0. The distance to be 

computed is 

L, = length AB = length FG = length FE~ length GE 
d 

L d = length FE - ift (U,18th HE) 

; If the level length is longer than the distance TES 

traveled by the submarine during a complete weapon search 

turn, the critical depth is at a level above L, but below 

d 

L/ . , ». Then the level defined by TES must be found and 
\d“l ; 

designated and the original level dropped from con= 

sideration until later. If L is shorter than TES then it is 

above the critical depth and its parameters must be computed 

to find qSV , , the volume searched in the d th level on the u th 
odu' 

search turn, (q = A,B). 
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TJ 




Figure 21: A sample volume portion for computation of L , . The points 

d 

(TJ^ d) are the radial endpoints of the corresponding 



planes. The level is under consideration and 



L 3> 1s 


to 


be computed. 


The 


bc 3 - 


SC 3 


= length 


FE 


BS 3 


ss 3 


= length 


BG 


ss 3 


ss 3 


= length 


HI 


SC 3 " 


sc 3 


- length 


HE 
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Note 8: See Figure 22. Subroutines RJ8, 9/ and 10 compute the distances 

and angles on the right portion of the figure which have not 

already been computed « Then subroutine RJ11 computes all of 

the special volumes in level After RJ11 the level index 

as well as the index r or s, is increased by L If index s 

is changed it means that the plane which determined the level 

was above plane TJ. Thus d and s must be increased together 

so that the difference is constant. It is their difference 

that is used for the comparison described by note 5° 

Note 9 l This subroutine computes the scanned volume in level 

(where d^D) on each search turn after the first . The first 

comparison is used to find out if the u th search turn is the 

last effective search turn in level L^ 0 If it is not* then 

subroutines RJ12, 15/ and 16 are used to compute the volumes 0 

The cross-section in this case looks like that of the right 

portion of Figure 22. If the u th search turn is the last of 

the effective search turns in level L,, then the normal cross® 

a 

section is cut at some point as in the left portion of Figure 22o 
In this case, the point where the cross®section is cut is 
found by one of the subsequent three comparisons and the cal- 
culation of scanned volume is modified accordingly 0 
Note 10: This corresponds to a vertical through the transducer which 

goes through a point in level such as point S in Figure 22 „ 
Note 11: This is used to compute the height of the critical depth 

above the search depth (length OR in Figure 22) and a correc® 
tion factor to be used to find the constants for level 
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Figure 22: Vertical cross-section through levels Lg and L ^ . 

The right portion depicts the volume scanned on the first 
search turn. The vertical OAB extends down to the trans- 
ducer . 

Volumes corresponding to depicted areas are 
Area Volume 



HID 

DEIJ 

EJK 

CFG = EJK 
ABCF ^ DEIJ 
HDEK 
ABCG 

The distances used to 
length NQ = TES 
length BD = 
length BE = B qo(d _ 1} 
length AH = B q ^ d 
length AK = B qod 
length HK = L d 



LqDV 

MqDV 



dl 



dl 



RqDV 



dl 



RqDV du ^ RqDV 



dl 



MqDV du t MqDV dl 



qsv 

qSV 



odl 



odu 

compute volumes in level L, are 

d 



length HI = b 



Id 



length JK -- length FG - b 



od 



length AB s h qd 
length AR ^ HL^ 
length BR = HL^ 
length DE - ^ 
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Note 12: This is the volume computation for FqDV^. This computation 

begins by turning the volume upside down so that level L, is 

U 

above level Then the computation is the same as that 

for (corresponding to area ABDH in Figure 22) required 

later . 

Note 13: See Figure 23. 

Note lh: See Figure 2b. 

length A I = D tan# 
area IJK = D 2 tan (6/2) 

volume AJKI = (l/3) (area IJK) (length Al) 

= (l/ 3 ) tan ( 6 / 2 ) tan# 



Note 15: 



Note 16: 



Note 17: 



Note 18: 



The volume computations (qMV^) for levels below the critical 
depth are begun at this point. 

The volume computations for the last level use B for the 

previous level and the angle j + d ■= r + p because there is 

no plane j on which to base the calculations. 

The calculations below search depth are identical to those 

above after the H. below are substituted for the H. above. 

J J 

The final section computes the scanned volume on the u th 
search turn separately for above and below search depth. 
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Figure 23: Volume computation of a triangular 









EB 



= A. Length BC = A tan (6/ 2). Length ED ■ h^ 
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•H 
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The altitude lengths are AI and AH 



APPENDIX C5 



FIRST TURN SCANNED VOLUME: i TH PING 




Note 1: p is now the angle between R q and the radial 

(R jk, 
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APPENDIX C6 



SCANNED VOLUMES SUBSEQUENT TO FIRST TURN: i TH PING 
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(?) |T Id -l^ijlEsTCT 



Note 7 



> 


( Note 8 




1 

1 



0 



(CH) 2 tan(&/2) [h qd - (B qI(d-1 j - (u-l)TES - CH)tan0 qId j 
- (l/3)tan(&/2)W> qia [(CH) 3 - (B qJfl - (u-l)TES)^] 



+ ASV— ASV 




(?) — ( 1 d d ) • • (?) 



Note 9 




|B_w, r (u-2)TES l CH 



Note 10 



n> 




t 


L B qIa - (u-2)t50h] 


i<J 


(Note 11 


1 ^ 







* ST 1 ST oa„^ S7 l O 



(l/3)tan(&/2)tan£ qId [(CH) 3 - (B H - (u-2)TES) 3 ] 

- qM? du + ASV -5* ASV 










'Note 11 




l> 





□ 



(CH) 2 tan(6/2)h d - qM7 du * ASV ASV 



|B qId - (u-l)TES : CH 
”> 



,< 


'Note 12 


I' 

1 





0 



(CH) 2 tan(b/2)[h qd - (B ql(a _ 1) - (u-l)TES - CH)tar>0 qld ] 

- (l/3)t«n(6/2)tan^ qla [(OH) 3 - (B gH - (u-l)TES) 3 ] * tSV-3-iSV 

3=5) 



©-C 



1 + d -5»- d 



L| ASV ^qSV lu ) ( 1 + u — >• u 




. f 



Note 14 




j' Note 13 



q I A>=— l 






<B^q 



D -^.d. 



D« -5»-d, 



0 

0 
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Note 1 
Note 2 

Note 3 



Note 4 
Note 5 



Note 6 
Note 7 
Note 8 



: p is now the angle between R q and the bisector of the ith ping* 

: In this case the target volume cylinder is large enough such 

that the scanned volumes after the first search turn are not 
effected * 

i In this case the target volume cylinder is small enough such 

that some of the scanned volumes after the first search turn 

must be modified * This is accomplished above and below the 
search depth and for each search turn by summing the modified or 
unmodified volumes over the levels to obtain qSV^* First 
each search turn above search depth is considered* then each 
search turn below * 

In this case the original qSV Q ^ is not changed <, 

: In this case the original LqDV^ and MqDV^ are not changed 

but RqDV^^ is modified and must be recalculated 0 The 
original cross-section appears in Figure 25 • Using formulas 
of subroutines RJ13 and RJl4 the new volume P is P ~ volume 
FGED + volume ADEC - volume ABGF * 

: In this case the original RqDV^ is zero* the original LqDV^ 

is unchanged and the modified MqDV^ is computed from the 
formula for subroutine RJ13« 

: In this case the original MqDV^ and R q DV du are zero and the 

modified LqDV^ Is computed from the formulas for subroutines 

RJ13 and 14 as the difference of two volumes* 

In this case the contribution of level L, on the u th search 

d 

turn is zero* 
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Figure 25: Computation of modified volume RqDV^ » Length FG is 

B q0(d-1) " ( u “ 1 ^ TES * Area BGH represents the original 
volume RqDV du „ Line AF is part of the vertical through 
the weapon. 
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Note 

Note 

Note 

Note 

Note 

Note 

Note 

Note 



9: The contributions of levels above the critical depth have 

now been determined* 

10s The contribution of level is not changed* 

11 s The missed volume on the previous search turn has been reduced 
leaving less to be picked up on the unchanged second turn* 

The formulas are those of subroutines RJ13 and RJl4* 

12s The missed volume in both the current and the previous 

search has been reduced* The formula is identical to that 
discussed in Note 7« 

13s The scanned volumes indexed by one set of q# u^ and i are 

now computed* The next step is to change u and compute the 
next set of volumes* 

lht Scanned volumes on the i th ping of all search turns above 
search depth are now computed* The next step is to compute 
the volumes below search depth* 

15: The increments of scanned volume for ping i ^ are now 

computed* The next step is to compute SV^* 

16s Having found SV iu> the index is changed to find 

SV,, , x . 

(i + l)u 
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APPENDIX C7 



TARGET VOLUME: i TH PING 
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Note 1: The problem is to compute areas as depicted horizontally and 

vertically lined in Figure l6 0 Using Figure 16 S with origin 

at point C and x axis toward point C, the equation of the 

target volume cylinder circumference on ping i is 

(x + R ) 2 + y 2 = R. 2 
x o' J 1 

Changing to polar coordinates (r^ p) with the same origin 

and measuring angles positive counterclockwise from the 

x axis, the equation becomes 

r 2 + (2R cosB ) • r=(R 2 =R 2 ) - o 

Solving for r, 

r = R. ^ 1 - (R /r. ) 2 sin 2 p =R cosP 
i 'o' i 7 o 

To find the area for ping i* r(r d p ) is integrated from 
p' + (i - 2)6 to p* + (i - 1)6 . 



Thus 

area = 

Let p* + (1-1)6 = B 
p» + (1-2)6 = A 



p' + (1-1)6 . 

T dp 

p*+ (1-2)6 



Then 
•B 
A 



f r 2 dp = f [(R* - R 0 2 ) ♦ 2 R 2 C0S 2 p 

-'a A L ° ° 5— 

- 2 Rq^V 1 - ( R o^i.) 8ln<t P C0S P 

J r B 

( R i 2 - R 0 2 ) dp = (R ± 2 - R o 2 )(B - A) 

A 



dp 
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2 R 2 cos 2 p dp = R 2 
o r r o 



p + 



(1/2) sin2p 



B 

A 



Let x - sinp 
Then dx = cosp dp 



-2R Q R i J" ^1 - (^Ai) 2 sin 2 p cosp dp 

= - 2 R^ f B Vl - (\A ± ) 2 7 dx 

•< A 



Let sinp = (R q /R^) x 
T hen cosp dp = (R q /R^) dx 



- 2 R^ j' \/ 1 - (R 0 /fe ± ) 2 x 2 dx = - 2 R ± 2 



2 „ 

cos p d p 



= - R £ 2 [p ♦ (l/2)sin2p 



B 

A 



Therefore * 




r 2 dp = (R q 2 - R ± 2 ) (sin2B - sin2A) / 2 
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APPEND IX C8 



ACQUISITION PROBABILITIES: I PINGS 



x , INPUT ... INEPT — _ 

1 ) ( APPENDIX C7 OUTPUT ) — ( Aw, W« , W» ) (1— » u p EI 




\/ 








sv lu PQi / rv ^ g ^ 


— ( i + i } — (i : m ^- 



L ( 1 U - g>U 




^ ~ 5>U g (M-l)u^ l// ^ ) ^ * g l(u*l)fc •“ j 




- ( 1 + u u V 



1 + 3 (1— ■ > U 



0-^g 



MU 






H g iu~ > g sl — ( 1 + s -^*0 — ( i ♦ i-»r 



1 + P ->U 





*-© 



(W'/Aw) - (W«/am) -5> R 


< 1 ^ r >-| 




1 



W' — =>■ w. 






[ (w lr “ w OSD^ PR V w lt 



I 



1 ->AP. 



[<f»s 



0H3>P A ( ''lr ) 



> £ 


-- x / 




N ip «^ ip 1 


— { i + i i ) — (i : ip 



1 - AP a -5> P A (« lr ) 





— i 


OUTPUT 


1 


( W’ “It? r = P - "•> R ) 



■0 
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APPEND K D 



PROBABILITY OF HIT 
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Note Is The problem constants are now determined. Next the 



probability P is determined for 
Hjk 



j = 1 > ° ° ° > ^i, > ^ = 1 , 



K. 
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